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INFLUENCE OF LARGE AMPLITUDES 
ON FLEXURAL MOTIONS OF ELASTIC PLATES 
By George Herrmann 

SUMMARY 


Starting with the fundamental equations of the general three- 
dimensional nonlinear theory of elasticity for the case of small elonga- 
tions and shears hut moderately large rotations , a set of plate equations 
of motion is derived for an isotropic material obeying linear Hooke's 
law. Certain assumptions as to plate displacements are made at the 
outset. The resulting nonlinear plate theory of motion, valid for large 
deflections, is discussed in the light of the three-dimensional theory 
and other nonlinear plate theories, in particular, the static Von Kannan 
plate theory. 

In order to investigate the influence of large' rotations in a dynamic 
problem, the nonlinear equations are solved for the case of propagation 
of straight-crested waves and the wave velocities are computed for various 
values of the parameters involved. Certain similarities to free vibra- 
tions of a mass on a nonlinear spring, governed by Duffing 's equation, 
are established. 

In addition, a variant of plate equations is derived which permits 
tracing bank the origin of certain terms appearing in other nonlinear 
plate theories, thus revealing more clearly the under lying assumptions 
of those theories. 


INTRODUCTION 


A variety of important problems of structural strength and stability 
of plates, arising in modern aircraft and missile construction, cannot be 
adequately analyzed on the basis of the classical theory because the plate 
deflections experienced are not small in comparison with the plate thick- 
ness. Several theories have been suggested to take into account the 
influence of large displacements and large rotations, the most widely 
known being the static theory of Von Karman. It is a common character- 
istic of these theories that the equations of equilibrium were obtained 
by considering the free-body diagram of a deformed plate element and by 
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equating to zero all the forces and moments. This procedure does not 
guarantee a rigorous consistency relative to the order of the terms 
retained, does not permit a clear identification of the order of the 
terms involved, and precludes the possibility of occurrence of plate 
stresses typical of a nonlinear theory which vanish in the linearized 
case. Moreover, since the relationship between) for example, the 
Von Karman theory and the general three-dimensional nonlinear theory of 
elasticity has "been established only as far as the strain and displace- 
ment components are concerned, a more complete study of this relationship 
is needed. In view of these shortcomings an attempt has been made to 
derive a large-deflection plate theory of motion, starting with the gen- 
eral equations of the three-dimensional nonlinear theory of elasticity. 

The idea of deriving approxl mate theories of plates and other 
bodies, one or two of whose dimensions are small as compared with the 
third, from the general three-dimensional equations has been, in recent 
times, repeatedly applied with success with in the class of linear theories, 
as for example, in references 1 and 2. It is demonstrated herein that 
the same procedure, offering analogous advantages, may be used in the 
case of a nonlinear theory. 

In a first section of the present report, the main equations of the 
general theory are written down in the formulation given by Biot in 
references 3 to 5- These are followed by plate equations deduced from 
the general equations by assuming at the start the same displacements 
as in the Von Karman plate theory and performing appropriate integrations 
in the expressions to obtain the potential and kinetic energies used in 
the formulation of Hamilton's principle. 

The resulting stress equations of motion appear to be different 
from those obtained by Von Karman and by Love, in that inertia terms 
have been added. However, the equations presented herein are reducible 
to Von Karman* s equations in the absence of external moments, tangential 
forces, and ccmpressional inertia forces. The present theory is also 
reducible to the classical plate theory for the case where the rotations 
are taken to be small as compared with the strains and the rotatory 
inertia is neglected. 

All the nonlinear terms of the present theory appear to be of the 
same type as in the general three-dimensional theory and may be desig- 
nated, using Biot's terminology, as "curvature" terms and as "buoyancy" 
terms. A "torsion" term, occurring in the three-dimensional theory, is 
absent in the plate theory. The plate stresses entering the equations 
are defined automatically during the process. 

The plate stress-displacement relations are also derived by starting 
with the corresponding relations of the three-dimensional theory and ' 
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performing appropriate integrations an d omitting certain terms. They 
appear to be precisely the same as those given' in the literature which, 
however, were derived quite differently. 

The present theory consists, in essence, just as in the linear case, 
of the equations of motion, the stress-displacement relations, and bound- 
ary conditions. The role of the latter, however, is greatly reduced in 
the nonlinear case. While, in the linear' case, the boundary conditions 
are understood to be those quantities, in number and combinations , which 
have to be specified at the boundary in order to insure a unique solution, 
the analogous conditions presently ascertain merely the nonviolation of 
the principle of the conservation of energy. 

As an application of the equations of motion of plates with large 
deflections, the propagation of harmonic, straight-crested waves in an 
infinite plate has been studied. The inclusion, in the theory, of the 
effect of large rotations introduces a coupling between the classical 
compressional and flexural waves, and a redefinition of the normal plate 
stress permits exhibiting the nature of the coupling term more clearly. 

If, for this purpose, the rotation is considered independent of the 
transverse displacement, the operator matrix of the set of two equations 
Is shown to be symmetric, just as it should be in any linear theory. 

The two coupled nonlinear wave equations are solved simultaneously 
and certain similarities to Buffing’s equation and his method of solution 
are noted in the process. The biquadratic velocity equation describes 
the propagation of two modes, one essentially compressional, the other 
flexural. It contains the effects of transverse inertia, rotatory inertia, 
longitudinal inertia, flexural stiffhess, and compressional stiffness. 

The influence of each of these effects on the wave velocity Is discussed 
in detail. 

Within the range of va l idity of the present theory, that is, for 
flexural wave velocities which are small as compared with the velocity 
of compressional waves, the influence of rotatory and longitudinal inertia 
is shown to be negligible. It Is remarkable that in this case the two 
coupled, nonlinear equations can be solved exactly. For a given wave 
length, the amplitude has the same effect on the wave length as the plate 
thickness divided by the square root of three. 

In order to gain a deeper insight into the assumptions underlying ■ 
the plate equations derived by Love (ref. 6 ) and the bar equations 
derived by Eringen (ref. 7) .> a variant of plate equations was derived, 
in the light of the three-dimensional theory, in the last section of 
this report. The same plate displacements were assumed, but more com- 
plete expressions for the components of strain were employed. 

Shortly after the .manuscript of this report had been completed, 
Reissner, in reference 8 , noted that Kirchhoff (ref. 9 ) had been the 
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first to analyze motions of plates with large deflections. The present 
equations coincide, in essence, with those given by KLrchhoff, who, how- 
ever, did not make an attempt to solve any special problems. The first 
to consider equilibrium of plates with large deflections appears to be 
Clebsch (ref. 10 ). Both Clebsch and KLrchhoff used different procedures 
and arguments in deriving their theories, as compared with the ones pres- 
ently employed. A detailed analysis of these differences would require 
a separate study and is beyond the scope of the present report. A crit- 
ical review of the work of Clebsch and KLrchhoff is found in reference 11 . 

Th i s investigation was carried out at Columbia University, the major 
part of it being under the sponsorship and with the financial assistance 
of the National Advisojy Committee for Aeronautics. The initial phase 
of the work reported was sponsored by the Office of Naval Research. 


SYMBOIS 


A,B wave amplitudes 

c phase velocity 

Cp phase velocity defined by equation ( 46 ) 


D plate stiffness 

E Young ' s modulus 


e xx-> e yz;°tc 
e yy} e zx.> a) y > 
e zz^ e xy^ cu z 


quantities defined by equations (3) 


external plate forces and moments defined by 
equations (18) 


f x ,fy,f z components of traction (external force, per unit original 

area acting on plate surface) in x-, y-, and z-directions, 
respectively 

G shear modulus; G = p 

h plate thickness 


m x ,xfty,q 

^Mn^N^Mns* 

Qn 
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K 

L 

^ 1 *^ 32 * 


kinetic potential 
•wave length 

plate stresses defined "by equations (68) 


^1^2->^3^12 I plate stresses defined by equations (l2) 

M1.M2.Mp2 f 


n.s 


Ql*,Q2* 

S 

So 

T 

t 

U 

u.v.w 


u.v.w 


Uo^Vo.Wo 


Vo 


W 

We 


x.y.z 


7 


coordinates perpendicular to and along the plate boundary, 
respectively 

plate stresses defined by equations ( 67 ) 
surface of plane face of plate 

surface of plane face of plate before deformation 

kinetic energy- 

time 

potential energy 

displacement components in x-. y-. and z -directions . 
respectively 

approximate displacement components of plate theory, which 
are functions of all three space variables 

plate displacement components defined by equations ( 7 ). 
which are functions of two space variables only 

plate volume before deformation 

elastic strain energy 

work done by external forces 

original coordinates of particle (before deformation) 
wave number. 2rt/L 
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€ lk> € 22' e 33 

6 25^ e 31 ,e l 2 


components of strain defined "by equations (4) or ( 5 ) 


v 


5,Tb£ 

P 


Lame's constants of elasticity 
Poisson's ratio 

final coordinates of particle (after deformation) 
mass density 


T 11' T 22' T 33 1 
*&> r 31’ r 12 J 


components of stress defined "by equations (6) 


A prime indicates differentiation with respect to xj a bar indicates 
an approximation a dot indicates differentiation with respect to time. 


BASIC EQUATIONS OF THREE-DIMENSIONAL NONLINEAR THEORY OF ELASTICITY 


The basic equations of the three -dimensional nonlinear theory of 
elasticity written down in this section are taken from Biot's work in 
the field (refs. 3 to 5)* A more extensive account of nonlinear theories 
of elasticity may be found in references 12 and 13 . 


Deformation 

Let the original coord in ates X, y, and z of a particle attached 
to the material before deformation become 


| = x + u 
n = y + T 
£ = z + w 


y a> 


after deformation. The Cartesian rectangular coordinate system Itself 
is considered to be fixed dinring this operation. The deformation of an 
infinitesimal region surro unding this particle is given by the linear 
transformation 
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cLg 

dt] 



dx + <2y 

qy 


Bu 

Bz 


-£ te+ ( 1 + £)* + £ 

= 8vr dx + ^ dy+ / L + BA 

Bx qy \ ozy 


dz 


dz 


dz 


Using the notation 


e xx 

II 

®yy 

II 


_ Bw 

e zz 

Bz 

e yz 

_ l/Bw 
2 l$y 


_ l/Bu 

e zx 

2\Bz 

e xy 

H|CM 

II 

Oix » 

l(<k 

2 W 

°V = 

l/Bu 

2\Bz 

(Sz = 

l/Bv 


^ ( 2 ) 


^ (5) 


and assuming these nine quantities to he small,, the components of the 
strain tensor, including terms of first and second order, may he defined 
as 
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e ll 

e 22 

6 33 

e 23 

6 31 

e 12 


exx + exy&z - ezxoy + + ay 2 ) 

e yy + e yz ai x " e xj^z + K^ 2 + “z 2 ) 

e zz + e zx^y “ e yz a5 x + + ^ ) 

^z + 2 ^ ( Szz ” e yy) + 2 a V e xy ~ \ CD z e zx - 2 ^V^z 

e zx + 2 °V ( e xx “ e zz) + 2 “ 2 CD x e xy “ 2 “z^bt 

e xy + \ “z ( e yy “ e xx) + | aye zx - J aye yz - \ c^ay 


L (4) 


These strain components are referred to a local coordinate system (l, 2 , 
and 3) originally parallel to the x-, y-, and z-directions, respectively, 
and undergoing the same rotation as the material at that particular point. 
This rotation is given, in first approximation, hy a vector with compo- 
nents ay, ay , and ay, which will he adequate if these components are 
small as compared with unity. 


The expressions (eqs. GO) for the strain components, which are 
valid for elongations and shears which are small as compared with unity, 
contain linear and square terms in the components of rotation. The 
linear terms in the components of rotation are, however, third-order 
terms as may he readily seen hy expressing the quantities e.y in terms 

of the components of strain and hy remembering that all quantities given 
in equations (3) are assumed to he of the first order (ref. 13/ p. 52 ). 


Thus, the expressions for the strain components simplify to 

Gil = + i(ay 2 + ay 2 ) 


e 22 


+ |K 2 • 

e 33 

- e zz 

+ K"/ 2 ’ 

e 23 

= Syz 

- i ayy 

e 31 

= e zx 


e 12 

= Gxy 

- \ 


> (5) 
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Even though the rotations and strains are small as compared with unity, 
they do not necessarily have the same magnitude . For this reason the 
squares of the rotation are retained. 


Strain Energy 

The variation of the total strain energy in an original volume V 0 
is given hy 

6W "JIL ( T H 5e H + T 22 5e 22 + T 33 5e 33 + 2t 23 Be 23 + 

v o 

21-31 Se 51 + 2 rrj 2 Be-J^jdx dy dz (6) 

where the represent the components of a symmetric stress tensor 

taken per unit initial area before deformation with respect to locally 
rotated axes (l, 2, and 3) and as a function of the fixed original 
coordinates (x, y, and z). 

The equations of equilibrium may be derived from the expression for 
the strain energy by the same process as followed in the linear theory. 

The components of strain are expressed in terms of displacement deriva- 
tives and the principle of virtual work, the static analogue of Hamilton's 
principle, is applied, followed by a partial Integration. 


PLATE EQUATIONS 


Consider a plane plate of constant thickness h referred to an 
xyz Cartesian coordinate system, the xy-plane being the middle plane of 
the plate and the z-axis being directed normal to that plane. The 
approximate displacements u, v, and w characterizing the plate theory 
with large deflections are taken as 


5 = u o(x,y) 


▼ = v Q (x,y) 


dw c 

Sx 

dwg 


} (7) 


= vb(x,y) 
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These displacements clearly consist of two parts , namely, the transverse 
displacement w Q , which characterizes the classical plate theory, and the 
displacements Uq and v Q , which serve to describe the displacements in 
the plane of the plate. The displacements given hy 'equations (7) are 
those of the Von Karmen plate theory (ref. 1 ^) as was pointed out by Biot 
(ref. 5). 


Using equations (7), the following approximations to the quantities 
given in equations (3) are obtained: 


e xx 


_ 8uq 3^w ( 


cbc 


- z 


_ ^o 


®yy 5 y 


- z 


8x 

V 


'ZZ 


= 0 


eyz ~ 0 
e xz = ® 


e = 

/-v 


■*y ~ 2 \sr ^r 


8 2 w g N 

22 SrWy 


- _ dw 0 

^ 8y 

dw 0 

^ 2\ax 


Sy / 


(8) 


Plate deformations may produce only large rotations ajj. and ay 
but not 2)2 . Thus, in calculating the components of strain from equa- 
tions (5) with the aid of equations (8), second-order terms containing 
cUj are neglected. 


The approximate components of strain in the plate are then 
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z 




2 



» = lfoo + 5^0 _ Pz d Vp ,^ 0^0 

■*-2 2\^x Sy Sx Sy Sy Sx , 


^ ( 9 ) 


Strain Energy 


Inserting the approximate expressions for the strain components 
given in equations ( 9 ) into the expression for variation of the strain 
energy (eq. (6)) gives the following approximate expression for the 
variation of the strain energy 5 W: 



( 10 ) 
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Since the z-depen&ence of all the components of strain has been made 
explicit by equations (7), an integration through the thickness h of 

the plate from z = - to z = ^ can be performed immediately with 
the result: 



the integration being extended over the original surface S D of the 
plate before deformation. In this process one is led automatically to 
the following definitions of plate stresses: 
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Noting that 5| Suq, 


=A 

^8x / 8x 


^ w o\ _ 8 2 


6w 0 , and so forth, it is 


v8x / 8x"" 

possible, in expression (ll), to bring out the factors 5u 0 , 6v 0 , and 

S'W'o by partial integration obtaining 


SW 


-IL 


8N n 


8 2 M 1 


8x 


Su o ^ N 1 t- 2 5v 0 - 


8x 


AL * 

&c\ 1 8x 


8n 2 

$y 


5v^ - 


S 2 ^ 

Sy 2 


8w_ 


8v ° • M 1 * 2 • 


A 

8x 


("5 ^) Bw ° - ^("3 ^®»o - 


dN n p 8 n- 

1? 5v ° 


6v„ - 6u 0 - A(n lo ^ )6Wo . 


$y 


feP 2 8y 


^( h 12 s 2 ) 8 ^ - 2 

/ 


dx dy + 


Nji 5u n + Njjg 5u s + N^ ^-2. 6w 0 + Nn K ^-—2 5w 0 - 


8n 


on \ 8s 8s / 


8s 


ds 


Herein, use has been made of 


-(j> Mns ^7 6w 0 ds = (j) 6w 0 ds 


(13) 


(14) 


Work of External Forces 

■4 e 4 fx^ fy^ snd f z denote the x-, y-, and z-components, respec- 
tively, of the external force (traction) per unit original area acting 
at the boundary. The (virtual) work done by these external forces is 

rv/\ 

(15) 


8W_ 


-IL <* 5u + fy 5v + f z 6w)dS 


the integral being extended over the whole surface of the body- 
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In the case of a plate, the total surface consists of the two plane 
faces z = and the cylindrical surfaces hounding the plate. The 
expression for the work is therefore 


■ [//( 


fy 5 u + f v SV + fn 


~ 1 z=h 
5 w)dx dy 

— 7.=- 


=h/2 


// J (f n 5 ^ + f s 5 u s + f z Sw + f n 5 w +' f s bvj dz ds (l6) 

The line integral is to he taken around all (external and internal) 
cylindrical boundaries of the plate. Subscripts n and s indicate 
components referred to coordinates n and s measured normal to and 
along the boundary. 

Substituting for the displacements u, v, and w their approxima- 
tions u, v, and w from equations (7), the work of external forces 
acting on the plate is given by 

_ ( PP \ r /8w„\l / . dwj\ ] _ \ 


" [jj ' fx ^° " Z 8 (sf)_ 

p r h / 2 T (< 

(b J f n 6uon - z 61 - 

J J -h/2 L 


f n ^ + *s |^ ) 8w 0 l dz ds 
on os / 


+ f y Bv 0 - z' — + f z 8w 0 dx dy 


B=-h/2 


+ fs 5 uqs - z 8 


+ fz 8w 0 + 
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Employing the notation 


/ s h / 2 

** = ( fx )-i/2 


h/2 

-h/2 

h/2 


Etc “ '(zfx) 

N n = / ; f n 

J -h /2 

• ./■'■ ,, 

v. /o 


dz 


= / fnZ dz 

% - (%) 

= ^. h /2 

> h/2 


-h/2 
h/2 

) 

-h/2 
.h/2 


Bna* -/J 


On 


■/: 


■h/2 
h/2 
■h/2 


f s dz 


f z dz 


^ "X 

expression (17) may he rewritten as 


1 = 

h/2 
h/2 


f 0 z dz 


(18) 


8 tL 


F x 8uo - mx 81 + Fy 8y 0 - ™y 8l ^— J + q. 5w o 


dx dy + 


■ -IT 

jf [%* 5uo n - Mn* 5^J + Nns* 8uo s - Mqs* + Qq 8v 0 


N n * ^ 5w o + %s* ^ 5w o 


ds 


( 19 ) 



16 


NA.CA TN 3578 


Noting that 


If m x 8^) + ry 8^) 8x dy 8v 0 ds - JJ + |^6v 0 dx dy 


the work of external forces may "be expressed as 


5W C 


F x 5u 0 + F y 8v 0 


+ L + ^ + ^l) 


5v r 


sr ^ 5T/ 0 


-1 


^ Nn* Su^ + Nns* 5u s - Mq* b(^j + ^ 


dx dy + 


8 Sw n - 


m n 6w 0 + Sw 0 + %* ^2. Sw Q + N^* ^ 8w 0 

on os 


ds 


(20) 


(21) 


Kinetic Energy 

The kinetic energy T of a body occupying before deformation a 
volume V 0 and having a mass density p is given by 

T = J JJf p(u 2 + v 2 + v 2 ) dx dy dz (22) 


The dot indicates differentiation with respect to time t. Using u, v, 
and i? in the above expression, an approximate form of the kinetic 
energy T in a plate is obtained as 


T 


-iJIL 


8w, 


,v \ 2 


8w 0 ' 


. 2 


• '-'«o l . l • '-'wo » . 

r si + r • 2 *rl + v ° 


dx dy dz (23) 


which can be integrated through the thickness h of the plate yielding, 
under the simplifying assumption that p is independent of z. 


T 


- 1 /To hu 2 + h 5 /^ o\ 2 y. 2 h 5 /^ 2 . x. 2 

- g JJ P + hvo + + 


dx dy (24) 
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The variation of the kinetic energy 8T is found to he 


8T 


■ If 


phiio 8u 0 + pht 0 8v 0 - 


3 


w„ 


Ph- 

12 2k 2 


5*,-, - 


8w 0 + phw 0 8w 0 ] dx dy + 

12 <*r 


/ 4 s 25 '” 48 


(25a) 


By integrating partially the above expression with respect to time and 
setting, as usual, the variations at the beginning and end of the time 
interval equal to zero, the following equation is obtained: 


5T 




huo 8 uq - 


hw Q 8w 0 dx dy 


J 12 Sn 


8w n ds 


(25b) 


Equations of Motion 

Hamilton's principle is now applied to derive the equations of 
motion. The principle states that for an arbitrary time interval t-j_ - t 2 


8 



Kdt = 0 


(26) 


where K is the kinetic potential} in general, 

K = T - U (27) 

U being the potential energy. In the present case of the plate, the • 
potential energy consists of two parts, namely, the approximate internal 
strain energy W and the approximate potential energy of the external 
forces -W e . Hence, 


K = T - W + W e 


(28) 


and 

8K = ST - 8W + 8W e (29) 

Using equations (25), (l3).> and (21) for 8T, 8W, and 8W e , respectively, 

the following expression for the variation of the kinetic potential is 
obtained: 


/ 
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SK = JJ P f-Mio Suq - bv Q 5v 0 + ^ Bw o + &w o “ 


h^o SVq) dxdy - a* +JJ 


V 8 2 Mq_ 

— 6llo + — ;r— &Wn + 


8n- 

6uq + — o~ Bw o 
^ 8 k 2 


l( H1 1 2 ) 5W ° + ? 6T ° + 5w ° + i( Ha &) 
I 2 ) ■* ♦ #(% £*) 


Bw 0 + 


&v ° + 5v ° + “§? 8U ° + 


^xV 12 8y 

§ 


bw. 


^-(w.o =-a] 5v 0 + |-f Nl2 — 


2 

Bw ° + 2 ST* Sw ° 


dx dy - 


N n 8u n + W ns 8u s + N n ^ Bw 0 + Hqs ^ Bw 0 - Mq 


8n 


F x 8uo + Fy 6v 0 + 


( 2 l ea + S a ) 6 wo ] 4 s+ i 5 r 

+ lx^ + 8w °| ^ ^ + ff ^n* 8u n + N ns* Bu s + M n* B (^) + 


Qn - %i 


* 8w n „ * 8w 0 

r- 2 - + 13ns ^ 
on os 


+ ^ns_ + it* ^ + Hns 

OS 


&W r 


ds = 0 


(50) 


This variation of the kinetic potential 5K must he equal to zero 
for arbitrary values of 5 uq , &Vq^ and so forth. The three stress equa- 
tions of motion are obtained from the integrand of the double integral in 
equation (30)^ obtained by equating to zero the coefficients of 8 Uq^ 

Sv 0 , and 5w OJ are 
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The integrands of the line integrals serve to identify the Btarred 
with the corresponding unstarred quantities (such as N n * = N n ) and permit 
the establishment of a relationship between the transverse shear force 
and other plate stresses as follows: 


On = 


8M n ^ns 
8n 8s 


+ m n + 


ph 3 

12 8n 


( 32 ) 


In addition, the integrand of the last line integral may be employed 
to establish the boundary conditions. For a system of linear differential 
equations, appropriate initial and boundary conditions are those which are 
sufficient to assure a unique solution. The classical uniqueness theorem 
of the linear theory of elasticity, using energy considerations, is given 
by Neumann (ref. 6, p. 176) . its analogy for a linear plate theory, 
including the effects of rotatory inertia and shear, was established by 
Mindlin (ref. l). 

Similar arguments may be used in the present nonlinear case. Since 
the superposition principle does not hold, the appropriate boundary con- 
ditions will ascertain not a unique solution but merely the nonviolation 
of the principle of conservation of energy. These boundary conditions to 
be specified on the edges of the plate are determined by any of the com- 
binations which contains one member of each of the four products of the 
last line integral in equation (30). 
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The first two equations of the system (eqs. (31)) describe motions 
in the plane of the plate and possess the same form as the equations of 
the plane stress problem. The third equation describes the transverse 
motions of the plate. The first three terms are those of the classical 
plate theory. The following six terms represent the influence of large 
rotations and are nonlinear. Each of these six terms , for example, the 

-v / dw \ 

first, -2- N, — 21 } may be vritten as the sum of two terms; that is, 

8x \ 8x / 


cix 




8w q 

8x 8x 


+ N, 



(33) 


Inasmuch as dw 0 /cix represents a component of rotation, it is 
recognized that the nonlinear terms in the stress equation of motion of 
the present plate theory are precisely of the same nature as those in 
the three-dimensional non li near theory of elasticity. Following the 

8w 

terminology of Biot (ref. 4) terms of the type — = — ^ may be called 

ox ox 

buoyancy terms, because they arise from some kind of buoyancy due to the 
deformation of an element in its own stress field and depend on the stress 

8 2 w r 


gradient, while terms of the type Nq 


have their origin in change 


8x 


of direction of an element due to curvature and may be called curvature 
terms. 


The two terms containing the total transverse normal stress will 

be neglected in conformance with the type of motions to be described by 
the plate equations for which this stress is negligible, just as in the 
classical plate theory. The remaining three terms on the left-hand side 
of the third equation of motion (eqs. (30)) are the forcing terms, q 
being the net normal pressure, while m x and my are the external moments, 
which are customarily neglected in a plate theory. They are retained here 
merely for the sake of completeness. 


On the right-hand side of the equation are found the inertia terms, 
the first pfaw 0 being the transverse inertia term of the classical plate 

Ph 3 . Ph? 8 2 w q 

12 ^2 12 x 2 
ox oy 

which are discussed in reference 1. 


theory while 


are the rotatory inertia terms 


There exist several other possible ways to deduce plate equations (3l) 
from the general theory, two of which are outlined briefly as follows: 
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First, the three-dimensional equations co uld he written down using 
the components of strain given by expressions (5), which would yield, 
after substitution of plate displacements (eqs . (7)) and appropriate 
integration, the plate equations (31). 

A second alternative consists in following the same procedure which 
leads to equations (31) hut temporarily consider ing the rotations uy 
and ay as being independent of the transverse displacement w Q in 

perform in g the variation of the energy expressions . This would lead to 
five equations of motion, with five unknowns u Q , v 0 , ay, ay, and w D . 

The equations would contain, in addition to the plate stresses defined 
hy equations (12), the plate shears Qx and Expressing ay and 

ay subsequently as derivatives of w Q and eliminating the shears Gy 
and Qy from the equations of motion lead to equations (31), two inter- 
mediate relationships in the x- and y-directions being those analogous 
to equation (32). 

Comparing equations (3l) with the Von Karman plate equations as 
found, for example, in reference 15, it is noticed that the essential 
difference consists in the fact that the Von Ka rman equations' do not con- 
tain the buoyancy terms, but only the curvature terms, as the nonlinear 
terms of the stress equations. Thus, a mere addition of inertia nnfl 
forcing terms to the Von Karman equations will not result in equations (3l). 
However, in the absence of inertia terms, equations (31) reduce readily 
to the ones of Von Kanaan. 

Comparing equations (31) with those given by Love (p. 558 of ref. 6), 
the absence, again, of the buoyancy terms is noted and, in addition, the 
occurrence of nonlinear terms containing the transverse shear force. The 
discussion of the significance of these latter terms 1 b postponed to the 
last section of this report. 


Stress-Displacement Relations 

The stress equations of motion (eqs. (31)) represent three equations 
relating 10 unknown functions, namely, the 3 displacements and the 7 plate 
stresses. Therefore, just as in the three-dimensional theory, additional 
equations, here numbering seven, are needed to complete the problem. They 
should relate plate stresses to plate displacements - and are developed 
from the three-dimensional relations, as was done before for a lin ear 
plate theory by Mindlin (ref. l) . Assuming the material to be isotropic, 
the stress-strain relations expressing Hooke's law are 
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T 11 " ^ ( e ll + £ 22 + e 33) + ^11 

t 22 = ^ ( e ll + e 22 + 6 33) + ^22 

t 33 = A ( e ll + e 22 + e 33) + 2lie 33 

t 12 = 2 ^ e 12. 
t 23 = ^ e 23 
T 3l = ^ e 3l 


f (5*t) 


where the stresses and strains are those Introduced by expression (6), 
and A and p are Larne's constants of elasticity. The third equation 
of the set (eqs. (3^-)) is now solved for anh. substituted Into the 

first and second, resulting in 



Ev 


e ll + 

! - v 2 ^ A + 2p 

T 33 

6 22 + 

Ev . A 


1 - v 2 ‘ L± A + 2p 

t 33 


Performing the integrations indicated in equations (12) and neglecting 


/ h/2 nn/ii 

Tit dz and, in addition, / T 33 z <3- z results In 
•h/2 -h/2 


h/2 


-h/2 


h/2 


the following plate stress-displacement relations after substitution of 
the strain-displacement relations in equations ( 9 ): 


M-i = -D 


'8% 

V8x 2 


+ v 


w r 


8y c 


Mg = -I 


■w> 


w r 


V8y 


chc 


S 2 w 

= -D(l - v) 0 


8x 8y 

(equation continued on next page) 
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N i = 


Hr 


Eh 

^0 . ll*«o\ 2 . v ^r 0 

1 - v 2 

dx £\dx / dy 


~ . 2 

,Eh 

^0 lj OW 0 \ ( ^ ^Uq 

H 

1 

*10 

$y 2\dy } dx 


+ vi 




2\8x 


/dir. dv Q dw Q dw Q \ 
N-, 0 = (2a —2. + —2. + — 2. — ° 

12 \dy dx dx by ) 


(35) 


In relations (35) E is Young's modulus, G is the shear modulus, 
v is Poisson’s ratio, and D is the plate modulus expressed as 


D = 


Eh > 


_ Gh > 


.12(l - v 2 ) *>(1. - v) 


(3 6) 


The stress-displacement relations (eqs. (35)) coincide in every 
respect with those of the Von Karman plate theory (ref. ]Jf) R-nfl also 
with those of the more general theory of Love (ref. 6), which are, how- 
ever, derived in a different manner without making use of Hooke ' b law 
for a three-dimensional solid. 


STRAIGEEC-CRESTED WAVES 


As an application of the plate theory derived above, the problem 
of propagation of free waves in an infinite plate is now investigated. 
Ass umin g the straight-crested harmonic wave to be propagated in the 
x-direction, the equations of motion are simplified by letting v Q = 0 
and d/dy = 0. 


In the absence of external loads the one-dimensional equations are, 
from plate equations (3l), 


Hq' = phiio 

M l" + %-( N l w o') = P^o “ V 


(37) 


Primes i n dicate differentiation with respect to x. The corresponding 
stresB-displacement relations are, from equations (35), 
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so that the displacement equations of motion take the' form 


(38) 


Eh 


1 - v 


r(V + VoW) = Pkuo 


-Dvr 0 1V + 


1 - V C L_ 


11 • 1 11 11 / i\ 2 •• 11 , .. 

Uq V 0 + Uq w d + ^ v q (w 0 ) + Wq = phw, 


(39) 


It is noted that the first of the two displacement equations of motion 
(eqs. (39)) describes essentially the compressional motions in a plate, 
which are, however, coupled to the transverse motions because of the 
presence of a (nonlinear) term containing derivatives of transverse dis- 
placement w 0 . The second equation governs the transverse (flexural) 
motions, which appear to be coupled to compressional motions by the pres- 
ence of two terms containing derivatives of the longitudinal displace- 
ment Uq . 

Neither the form (eqs. (37)) of ft 16 stress equations of motion nor 
the form (eqs. (39)) of the displacement equations of motion permits 
clear recognition of the nature of the coupling effect contained in these 
coupling terms. A deeper insight into this coupling phenomenon in both 
sets of equations ( 37 ) and ( 39 ) may he gained, however, with the aid of 
a redefinition of the plate stress N-^. 


In view of the stress-displacement relation contained in the first 
equation of the set (38), N-|_ may be expressed as the sum of two compo- 
nents, namely, 

N i = «x - fro) 


where N x is the classical component defined as 


Nx = 


Eh dup 
1 - v 2 dx 


frl) 


and N zx is defined as 


^zx 


Eh 

2(l - v 2 ) 8x 


(42) 
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The rotation component ay appearing in equation (40) is considered, 
temporarily, as being independent of the transverse displacement w Q . 
The resolution of U-j_, defined hy equation (40), is convenient because 
it separates the terms containing Uq and w Q . In terms of the newly 
introduced plate stresses and N zx , the stress equations of motion 

(eqs. ( 37 )) W be written as 


Nx - j^(NzxPy) - phu D 


M l" - = P^o - ^ ^ 

cbc 


and the displacement equations of motion (eqs. (39) h, as 


Eh „ >■ 


Eh 


ctf. 


rr* v ax* sf; - phu ° 

iv 


-Dw, 


Eh 


Eh _d_/ ^o ) 

1 _ v 2 dx / 2(1 - v 2 ) 8 x 


b( 2 ^o\ 

“S — ) = P hw o 


> ( 43 ) 


cbc 


12 


$ .. 

W r 


ph' 11 


> m 


The coupling terms, the second term in each of the equations of the 
set (43) or (44), are thus clearly exhibited. As a consequence of a 
large rotation ay, the plate stress N x enters the equation of trans- 
verse motions, while the plate stress N zx , which depends on w Q only, 
enters the equation of compressional motions. 

The coupling operator in equations (44) is recognized to be the 
same in both equations, namely. 


1 


Eh 

- V 



. Eh 

1 - v 2 ^ 



( 45 ) 


thus making the operator matrix on the two dependent variables Uq and 
w 0 symmetric. Thus, it appears that the operator matrix is symmetric, 
in a generalized sense, even in the present nonlinear theory just as it 
should be in any linear theory, as was pointed out, for example, in 
references 1 6 and 17- Moreover, it is observed that the coup ling oper- 
ator contains both a curvature term and a buoyancy term. 
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In addition, it is noticed that the effect of a large rotation is 
apparent not only in the coupling terms, -which contain the rotation to 
the first power, hut also in a term occurring in the second equation, 
which contains the rotation to the second power. 

Equations (39) are now rewritten in a more convenient form. The 
first is solved for w 0 ’w 0 " and substituted into the second. If Cp 
denotes the phase velocity of compressions! waves in a plate, then 





and equations 


(39) may he written as 



(46) 


(47) 


The solution of this coupled set of two partial, .homogeneous, nonlinear, 
differential equations is sought in the form of harmonic waves 


Uo(x,t) = A sin 7 (x - ct) 
v 0 (x,t) = B cos 7 (x - ct) 


> (48) 


where A and B are the amplitudes, 7 and 7, the wave numbers, and 
c and c, the phase velocities. Suhstitution of this trial solution 
into the first of the equations of motion reveals that, for periodic 
motions, 

7 = 27 (49) 

c = c (50) 

8A7 ^1 - sLj = B 2 7 2 (51) 

The first term in equation (51) represents the influence of compressions! 
Btiffness, the second, the influence of longitudinal inertia, and the 
term on the right-hand side, the influence of transverse stiffness. A 
sketch of the deformed plate is given in figure 1. Suhstitution of equa- 
tions (48) into the second equation of the set (47) results in the 
following equation: 
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2 h^ . 

“ C P H ^ 


cos 


A + 2ky^j + c 2 {^k-p + + 7 2 j 

^ky^c 2 cos 37 (x - ct) =0 

In obtaining equation (52), use was made of the relation 


7 (x - ct) - 


(52) 


COS' 


„ _ 1 
a — 


cos 


3a + ^ 
4 


cos a 


(53) 


Equation (52) is similar to an intermediate form of the Duffing equation 
which governs, for example, motions of a single mass attached to a non- 
linear spring and -which is discussed in detail in reference 18. Discus- 
sion of nonlinear problems governed by partial differential equations 
may be found in references 19 to 21. 


Confining attention to periodic solutions, the term (eq. (52)) with 
cos 37 (x - ct) may be disregarded as compared with the term 
cos 7 (x - ct), provided that the amplitude of the former is small as 
compared with the amplitude of the latter. This restriction is discussed 
in detail below. The velocity equation is then 

4r(i + + »r) + 2 a 7 (510 

Cp 


The first term on the left-hand side of velocity equation (54) 
represents the influence of transverse inertia on the wave velocity, the 
second term, the influence of rotatory inertia, and the third term, the 
influence of longitudinal inertia. On the right-hand side, the first 
tern represents the influence of flexural stiffness and the second term, 
the influence of campressional stiffness. 


Eli m inating kj in equations (51) and (54), a single biquadratic 
equation 




+ 1 


+ + i?£ = 

12 4 


(55) 


governing the phase velocity c is obtained, which is seen to depend on 
two parameters, namely, hy and B7. Since 
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■where L is the -wave length, h 7 represents the ratio of plate thick- 
ly 

ness to wave length, while B 7 equals the maximum rotation oiy = - ^— 2 . 


This velocity, or frequency, equation represents a typical two-mode 
coupled system similar, for example, to the one of reference 2. While 
in reference 2 the coupling "between the two modes was established through 
Poisson's ratio and thus depended on material properties, the coupling 
is effected presently through a ki n ematical quantity, namely, by taking 
into account large rotations. A coupling due to gravity between the 
same two types of waves has been discussed in reference 3- The two modes 
may be uncoupled by omitting the terms in equation ( 55 ) containing B 7 , 
which yields 




c 2 h^7^ 
+ 7J~12- 
°v 



(5 6 ) 


Thus, one mode, which describes compressional motions in the plate, has 
the velocity 



(57) 


while the second, the flexural mode, has a velocity 



.2 2 

h 7 


1 + 


12 

.2 2 

h 7 

12 


(58) 


which has been discussed in reference 1. If the Influence of rotatory 
inertia is neglected, equation (58) simplifies to 



_ b 2 7 2 . 
12 


(59) 


Several graphs have been plotted to exhibit more clearly the coupling 
effect. Figure 2(a) shows the uncoupled and coupled velocity curves for 
a constant value of. the ratio of transverse amplitude B to the wave 
length L, which is a measure of the rotation and which has been chosen 
to be equal to 0.1 in the case where the rotatory inertia is neglected. 
Figure 2(b) gives the same curves, including, however, the effect of 
rotatory inertia. It is seen that the coupled curves approach the veloc- 
ity curves of uncoupled motions for large values of h/L but not for 
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h/L = 0. This deficiency is due, most probably, to the fact t hat the ' 
equations derived are not valid for the full range of the parameters 
involved. These inherent limitations are discussed subsequently. 

Noting that — = £ “ similar curves can be drawn for a constant 
L n L 

value of B/h, which may be taken as a measure for the deflection. Fig- 
ure 3(a) shows the two coupled velocity curves, drawn for 'S. = l^ with 
the effect of rotatory inertia being neglected. The coupled curves 
approach the uncoupled curves for y = 0 but not for large values of 

Xj 

h/L. Figure 3(t) gives the coupled and uncoupled velocity curves for 
the same constant value of B/h, but with the effect of rotatory inertia 
included, and a gai n a similar behavior of the coupled curves with respect 
to the uncoupled ones is noted. 

The range of validity of the velocity equation, however, is limi ted 
to specific ranges of the parameters involved in view of the simplifying 
assumptions underlying both the plate theory and the solution of the 
wave-propagation problem.- These limitations are now examined in detail. 

At the outset, in considering the possible deformations, the rota- 
tion was assumed to be not too large, so that it could be represented 
by a vector with components ay, and a^. This is admissible only 

if these components are small as compared with unity. 

In the case of waves, described in equations (48), propagated in a 
plate, the maximum value of the component of rotation uy is given by 

1^1 max ~ ^ 

which has to be small as compared with unity. Moreover, it was assumed 
that the components of strain are small as compared with unity, which 
implies that in the expressions for the str ains all the nonlinear terms 
are of third and higher order, except the squares and products of rota- 
tion, which are of the second order. Such a typical nonlinear term not 

vanishing in the one-dimensional case is i aye-^, which occurs in the 
expression for the strain e^]_ the set of equations (4). Another 

*i 0 

such term is — e^ , which occurs in the expression for the strain e-, - [ 
if the strains are not assumed to be small (ref. 4, eqs. (l-12)). 

Thus, the assumption of small strain requires the smallness, as 
compared with unity, of the coefficient e^ h In the case of the 
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present plate theory, the coefficient 


s 8u 0 

exx _ ^ is 


85 = ^ 

8x 8x 


- z 


>2 
O W,- 


and consists of a campressional part 


8x 

8x 


and a flexural part z 


the maximum of which is 


8 2 w r 


8x 


8 2 w, 


o 

2 * 


Since the coefficient 


8x 

exx has to he 


small throughout the thickness of the plate, each part must he small as 
compared with unity- In the case of waves, this Implies that hoth A7 

(compressional part) and — B7 2 (flexural part) should he small as com- 
pared with unity. 


The significance of the smallness of the flexural part is not quite 
clear yet, in view of the fact that B7 is small hy itself. However, 
it should he pointed out that, in the expression for the strain, the 


o o 

squares of rotation contain terms with B 7 which are retained, while 
terms of the type 1 uyexx which contain ^ B 2 7 2 are neglected. 


Thus, clearly, this implies that hy is small as compared with 
unity, and since 7 = 2rt/L, where L is the wave length, that h/L is 
small as compared with unity. Hence, the solution will he valid only 
for waves which are long as compared with plate thickness. 

To summarize the restrictions of the wave solution inherent in the 
present plate theory, the smallness of the components of strain implies 
the smallness of A 7 and I17, while the smallness of the rotations 
implies the smallness of B7, all compared with unity. 

Attention is drawn next to Limitations inherent in the neglect of 
the term with cos 37 (x - ct) as compared with the term with cos 7(x - 
in equation (52) . This neglect is justified only if the amplitude of the 
former is small as compared with the amplitude of the latter. The ratio 
of these amplitudes, which ought to he small, is calculated to he equal 
to 3 A 7. But this requirement coincides with one of the requirements 
imposed hy the smallness of strain discussed previously. Thus, the 
approximations introduced in the wave solution do not introduce any addi- 
tional limitations. 
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The solutions (51) and ( 5*0 are now examined in the light of the 
limitations just discussed. Solving equation (51) for A7 and obtaining 


A7 


= 


8 



it is noted that 1— — — should he small, as compared with unity, 

1 - c 2 /c p 2 

for a given value of B7. The curve of A7, as a function of c/cp, 
exhibits a resonance phenomenon of a type similar to that in the case 
of forced vibrations of a single mass. The phase velocity c should 
be either smaller or larger than the campressional wave velocity in a 
plate Cp. 

.2 2 

In equation ( 5*0 the terms with 7 - and 3A7 on the left-hand 

side may be omitted, in accordance with the limitations imposed, since 
these terms are added to unity. Thus, relation (3 & ) simplifies to 

4* = tg- + 2Ar ( 60 ) 


That is, in the range of its validity, the phase velocity depends only 
upon transverse inertia and the flexural and campressional stiffness, 
while campressional and rotatoiy inertias are negligible. But it has 
already been established that I17 and A7 are small as compared with 

unity j hence, from equation (60), is small as compared with unity- 

C P 

Thus, equation ( 5 l) simplifies to 

8A7 = B 2 7 2 (6l) 


Eliminating ky from equation (60} with the aid of equation (6l) gives 
the velocity equation 





+ B = 


=- + B 


(62) 


in which each term has to be small as compared with unity. The result 
is thus obtained that, within the range of its validity, for a given 
wave length, the velocity is affected in the same mann er by the deflec- 
tion B as by the plate thickness divided by \J 3 . 
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The velocity given by equation (62) is plotted in figure 4 for 
various values of the parameter B/L, which represents a measure of the 

rotation. Since ^ equation (62) may he rewritten as 

L h L 





( 63 ) 


where the parameter B/h may he taken as a measure of the deflection. 
The velocity given hy equation (63) is plotted in figure 5 for various 
values of the parameter B/h. 


In the treatment of nonlinear systems with one degree of freedom, 
it is customary to study the so-called response curves, that is, curves 
relating the amplitude to the frequency (ref. 18 ). Since, in the pres- 
ent case, the velocity may he said to play the role of the frequency, 
B/h was plotted as a function of c/cp for various values of h/L, 
which is now considered a parameter. These curves, given in figure 6, 
are qualitatively the same as the response curves of systems governed 
hy Duffing's equation (ref. 18 ) which demonstrates still further a cer- 
tain resemblance of the present nonlinear continuous system with infi- 
nitely many degrees of freedom to a nonlinear single-degree-of -freedom 
system. 


It is worthwhile to point out that, if the terms containing longi- 
tudinal inertia are neglected in equations (bj) by letting U 0 = 0 , 
equations ( 6 l) and (62) represent exact solutions of equations (V7). It 
should also he observed that the contribution of the buoyancy term in 
the second equation (^7) is contained only in this longitudinal inertia 
term. 


A VARIANT 1 OF PLATE EQUATIONS 


In this last section of the present report, a large-deflection plate 
theory is derived the equation for which includes the influence of terms 
containing the first power of the rotation in the expressions for the 
components of strain, that is, expressions ( 4 ) for the strains will be 
used and not expressions (5). Even though it has been recognized pre- 
viously that these terms are of higher order than the second, the devel- 
opment is worthwhile because it appears that better insight is gained 
in appraising other large-deflection plate and beam theories, particularly 
those derived by Love (ref. 6 , p. 558 ) and Eringen (ref. j). 

In view of the fact, however, that the assumptions as to displace- 
ments are still those given by equations ( 7 )> that is, e^ = e^ = e zz = 0 



NACA TN 3578 


33 


and. ooz is small, tlie influence of linear terms in cox and. cvy appears 
only in two components of strain, namely, e^-j_ and. eg which are now 


e 23 " 2 + 2 “V “ £ ° z ) 

e 31 = \ " \ ( e xy + 


> m 


or, vising the displacements given hy equations ( 7 ), 


2 

= - z 



- z 




1 8w 0 f ^Uq 

2 <ix\<dy 



_ 1 ^o ( 8w 0 
2 <dy \Sk 


- z 


-\2 
O W r 


Sx cty-y 


(65) 


Assuming, as was done before, that the derivatives of Uq and v Q with 

respect to the space variables are small as compared with unity and using 
the same procedure which led to the equations of motion ( 31 ), the following 
homogeneous equations of motion are obtained: 



In these equations, N-j_, Ng, N-|_g, M-^, Mg, and M-]_g are the plate 

stresses defined by the integrals in equations (l2). The transverse shear 
forces Q-]_*, Qg* are now defined by 
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%>* 



!> ( 67 ) 


and no longer by equation (32) as in the previous theory. Moments of 
transverse shear stresses M^]_ and are automatically defined as 



and do not occur in any linear plate theory. The static terms in equa- 
tions (66) could have also been obtained by a direct integration, through 
the thickness of the plate, of the corresponding three-dimensional equa- 
tions, which were given by Biot (ref. k, eqB. ( 5 A) )• 


Comparing equations (66) with those 
that the first two equations of each set 

that in reference 6 only shear curvature 

occur, while equations (66) contain both 
ancy terms, the latter being typified by 


given by Love, it is recognized 
are almost identical, except 

terms of the type Q-j* 


8x 

shear curvature and shear buoy- 
dQi* 8w 0 _ 
cbc etc 


the term 


Thus, the origin of the nonlinear terms in the first two stress 
equations of Love's theory is due to a retention of first-power terms of 
rotations in' the expressions for the components of strain. It has been 
shown previously in this report that they are of a higher order than the 
second and may thus be omitted. In reference 6 , it is merely stated 
that it will generally be sufficient to omit these terms. 

The last equation of the set ( 66 ) contains, in addition to the terms 
found in the corresponding equation of the set ( 3 l), a series of terms 
containing the moments M 23 and Mj-j_, which are missing in Love's equation. 
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In order to compare Eringen's beam equations (ref. 7 ) with those of 
the present study, the plate equations ( 66 ) are reduced to beam equations 
by letting v Q = 0 and d/dy = 0. They are 


ox ox \ ox 


£ * £(* &) 


^ fel 

Sx\ Bx Sx y 


= phtf 0 


Next, Eringen's equations are reduced to describe small strains and 
moderately large rotations by considering du/cbc to be small as compared 
with unity , letting the cosine of the angle of rotation equal unity, and 
assuming that the sine equals the angle itself. They are, in the nota- 
tion used herein. 


Sx dx 








* 8 ^ 


= Phw r 


Thus, it is seen that Eringen's term Q * ^ 1 ° is expressed in the present 


* Sx 2 


theory as 


^%1 


dx dx 


Columbia University, 

New York, N. Y. , October 3 > 1955* 
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(a) Classical theory, (l) flexural waves, (2) longitu d inal waves; coupled motions according 
to theory presented herein (b/L = 0.1), (2) lower branch, (5) upper branch. 




Figure 2.- Plot showing dependence of wave velocity on wave length and influence of large rotations. 
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longitudinal waves; coupled motions according to theory presented 
herein (B/h “1), Q) lower Branch, 0 upper branch. 


Figure Concluded. 
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